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We consider the Cooper-problem on a two-dimensional, square lattice with a uniform, perpendic- 
ular magnetic field. Only rational flux fractions are considered. An extended (real-space) Hubbard 
model including nearest and next nearest neighbor interactions is transformed to "fc-space" , or more 
precisely, to the space of eigenfunctions of Harper's equation, which constitute basis functions of 
the magnetic translation group for the lattice. A BCS-like truncation of the interaction term is per- 
formed. Expanding the interactions in the basis functions of the irreducible representations of the 
point group C^i, of the square lattice simplify calculations. The numerical results indicate enhanced 
binding compared to zero magnetic field, and thus re-entrant superconducting pairing at extreme 
magnetic fields, well beyond the point where the usual semi-classical treatment of the magnetic field 
breaks down. 

PACS numbers: 71.10.Fd, 74.20.Fg, 74.20.Rp 



I. INTRODUCTION 



The relationship between an external magnetic field 
and superconductivity is very important both from a the- 
oretical and a practical point of view. It is well known 
that superconductivity will disappear above a certain 
magnetic field due to orbital frustration within a semi- 
classical treatment of the problem, i.e. the pairing states 
(— fc, k) are not good quantum numbers. On the other 
hand, it was realized some years ago that there may 
be additional interesting effects as a magnetic field is 
increased to large values at low temperatures. Rasolt 
and TesanovicEl showed that re-entrant superconductiv- 
ity may appear for fields above Hc2 due to Landau level 
quantization. Because electrons in Landau levels will not 
be influenced by orbital frustration in a magnetic field, an 
increased magnetic field may in fact enhance supercon- 
ductivity. The re-entrance has no direct relationship with 
the low field superconductivity, and some materials may 
even be superconducting only in high magnetic fields! 
When Zeeman splitting was incorporated the re-entrant 
superconducting phase disappeared at high enough fields. 
We will investigate if a similar phenomenon also will hap- 
pen when a periodic lattice is incorporated, modifying 
the Landau level piclure of the single-particle spectrum 
to Hofstadter bands.B The infiuence of Zeeman splitting 
(i.e. effective Lande g-factor) on the Cooper binding en- 
ergy, will also be considered. 

In order to do this, we have considered the Cooper- 
problem on a two-dimensional, square lattice with a uni- 
form perpendicular magnetic field, i.e. we consider pair- 
ing of fermions in large magnetic fields at zero temper- 
ature. The same problem in zero field has been studied 
in Ref. ^. The model has on-site, nearest neighbor, and 
next-nearest neighbor interactions, where the values used 
for the interactions are input to the model. We do not 
have in mind any particular microscopic mechanism caus- 



ing these effective interactions, but assume that they can 
be attractive. Since the interaction potential has finite 
range, it is possible to expand the Fourier transformed 
version in a finite number of functions. These are basis 
functions of the irreducible representations of the point 
group C4iy of the square lattice. 

In the absence of interactions, the Schrodinger equa- 
tion for the problem is the well known Harper's equation, 
whose spectrum and eigenfunctions have been extensively 
studied, see e.g. Ref. 

The spectrum is the so called Hofstadter's butterfly. If 
the ratio of the flux per plaquette to the flux quantum 
is a rational number (^ — the energy bands split 
into q sub-bands. For irrational ratios the spectrum is a 
Cantor set.B Effects of electron correlations on the Hofs- 
tadter spectrum have been studied by several authors, see 
e.g. Ref. |l|,^9|j2^,|l|,|2|j2|, however they do not consider 
the possibility of Cooper pairing through an attractive 
effective interaction. 

For low magnetic fields the behavior of the system can 
be estimated by treating the field in a perturbative man- 
ner. Using the zero field correlation length one can find a 
value for Hc2 ■ The limitations of this approach is that it 
neglects the change in the wave-functions of the Cooper 
pairs due to the magnetic field. For very large magnetic 
fields, a more appropriate approach to the problem is to 
solve the single-particle spectrum first in the presence of 
the magnetic field, and then introduce the pairing in- 
teraction as a singular perturbation on this problem. If 
one is primarily interested in the behavior near Tc{H), 
one can exploit the fact that the order parameter/energy 
gap is small and develop linearized formulas as done in 
e.g. Ref. |1,||,|6[ 

We transform the extended Hubbard model (real- 
space) to the space of eigenfunctions of the Harper's 
equation which diagonalize the single electron part. A 
BCS-like truncation is then performed and the interac- 
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tions are expanded in the basis functions of the point 
group 6*4^ of the square lattice. From this Hamiltonian, 
an exact solution of the Cooper problem is in principle 
possible to obtain, however one encounters the problem 
that the associated matrix problem increases dramati- 
cally upon a reduction of the magnetic field. In fact, 
the matrix dimensions will be ~ x and only fields 
far above experimentally achievable ones can be handled 
in explicit numerical solutions at present .E2I We will con- 
sider extremely high fields where the flux per plaquette 
is a sizeable fraction of 0oi typically q < 8. 



II. THEORY 

A. The model 

We use an extended Hubbard model with nearest and 
next-nearest neighbor hopping matrix elements t and t' 
respectively, given by 
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Here U, V and W are effective interactions between elec- 
trons at the same site, nearest neighbor sites and next 
nearest neighbor sites, respectively. The noninteracting 
model including next nearest neighbor hopping has been 
studied in e.g. Ref. |2^. The Zeeman splitting is given by 
e{a) — —gsHBBa/2 where a — -l-(— )[=T ds is Lande 
(/-factor, and|-« is the chemical potential.E2l The Peierls 
phase factorsca for hopping from lattice site j to i is 
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where (jjQ = h/e. Employing a Landau-gauge A = Bxiy, 
introducing g — 2Trp/q — 2TrBa^ /(j>Q, and setting a = 1, 
the phase factor can be written as 
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g ^^^^^ n if \ri — Tjl = mex + ncy. 
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Using this, and writing the wave function as (we have 



chosen the normalization X^i^i Wti.v,iixi)\'^ = 1 where N 
is the number of lattice sites.) 



.. . 

where we have introduced k = [ji^v), the Schrodinger 
equation can be written as 

-e'^' {t + 2t' cos[g{m + 1/2) + t^)]} u^,,.,«(x„+i) 
-e"*'" {t -f 2t' cos[5(to - 1/2) + u)]} u^,^,£(a;„,_i) 
- {2tcos[gm, + v]}Uf^^^,i{xm) = f-ii.u,tUf^,„i{xm). (5) 

Here, x is taken modulo q due to the periodicity for 
m m + q, such that the functions u constitute the 
periodic part of Bloch functions on the magnetic lattice. 
The equation can be written as a g x q matrix. It is fur- 
thermore easily seen that the equation is periodic for v 
v+2tt. By introducing u^_^_f(a;m) = e*''^" ""-^ '"u^,y,^(a;m) 
it is easily shown that the equation for u^^i,^e{xm) is peri- 
odic for fj. ^ iJ, + 2TT/q. As a consequence of the fact that 
Xm - Xm mod q = nq, it is clcar that ^^.i,,e{x.m,ym) must 
have the same periodicity. The different eigenvalues are 
numbered with £ and each value corresponds to a differ- 
ent Harper- or Hofstadter-band. Using the completeness 
relation 1 = J2k. e\'^'^)('^^^\' ^^'^ creation/annihilation 
operators are seen to transform aS 
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where xi = xi mod q. Inserting Eq. ^ in Eq. ^ we get 

^0 = E ['^'^'^ + ^("^^ ~ 'i.Aa'^'^A'y (7) 
The two-particle terms in Eq. |l^ are all of the form 



i,j,{T,<7' 



By using Eq. || and a form of BCS-truncation of the in- 
teraction term (See appendix ^ for details) 



(T = — (T 

k' = — K, (— 7r/g, — tt) < < (tt/q, 7r), 
the two-particle terms can be transformed into 
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where 
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= E E^''^^.G.G'(^l,^2,A^)i?,,G,G'(^'i,4,'^')- 
G,G' Jj=l 

(11) 
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Here, we have introduced A = {U, V, W, V, W) and func- 
tions B given by 

B,j,G,G'{h,h,K) = u^^,,{G)u_^^i,{G')B,^{K), (12a) 



where 



and 



Bi{k) = y/N 

B2{i^) = Vn [cos(k2;) + cos(kj,)] 
i?3(K) = V2N cos(kj;) cos(kj,) 
Bi{K) = ^/n [cos(k2;) — cos(kj,)] 
Bc,[k) = V2N sin(Ka;) sin(Kj,), 
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(12c) 



Here G = {G^,Gy), where G^ = 0, 27r/g, ■ ■ ■ {q - l)2n/q 
is a reciprocal lattice vector for the magnetic lattice, and 



Gy^O. 



B. The Cooper- problem 

The two-particle Schrodinger equation is given by 

iHo + Hi)\l,2)=E\l,2). (13) 
Without interactions the Schrodinger equation is 
Ho\K/i,a;-Kj2,-cr)o 

(14) 

and it is possible to expand the two-particle states in the 
complete set | )o: 



eF)6'(e_K,,, 



^f) 



xaK,;^'l'«:^i><^;-«:^2,-o-)o, (15) 

where is the Fermi-level and 9{x < 0) = 0, 6{x > 0) = 
1. Inserting this in the Schrodinger equation Eq. |l^, mul- 
tiplying from the left by q{k,' , l'-^] — k,' , i'2\, using Eq. 
and introducing the shorthand notations x — {'^j^i,^ 
and ^ = {-q, G, G'} give 



(, X>XF 



(16) 



which implies 



-, X'>XF. (17) 



Here the notation x' > Xf means tK.>^t^^cr' > 
ep and e_K'.f^,-cr' > e_F- Inserting the expression for 
a^'.cr' into Eq. |l^, multiplying by B^i{x') and summing 
over x' > Xf gives 



(12b) 

where 



E ^I(x')i?«'(x'), 



and 



(19) 



(20) 



X>XF 

The final equation is 

E^«',«(i^,^'Me,.' =Ere,«^«,-'' (21) 

which has a (nontrivial) solution if and only if dct(D(i?) — 
r) = 0. This determines E and we can find the vector 
components ^4,0-', and thereby the Cooper wave func- 
tions in (k, ^)-space, a^'^f . The real space wave function 
is then 



(u|l)2)= E ^(eK,«i.<T - eF)6'(e-K,^2 -<T - ef) 

xat^- {i\K,h,a)o{j\-K, £2,-^)0, (22) 
^ ^ ^ 

or in mass-center and relative coordinates X — (xi + 
Xj)/2, p^ri-rf. 

(X,p|l,2)= E 0{€^j,,, - eF)6{e-^,i^,-a - ^f) 

K,,£i ,£2 ,cr 

X ai]-,'^u^MX + ^)u^^.,M - f )e"^-''. (23) 



III. RESULTS AND DISCUSSION 

In the previous section we have developed a framework 
for computations involving both nearest- and next near- 
est neighbor hopping and on-site, nearest neighbor and 
next nearest neighbor electron interactions. As a first 
step, we have studied explicit numerical solutions for the 
situation with only on-site attractive potential. This cor- 
responds to setting U < Q and F = = in the above 
formulas, i.e. the Hubbard model which in the absence 
of a magnetic field is expected to exhibit s-wave super- 
conductivity. The next nearest neighbor hopping term t' 
is also set to zero. 
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We have calculated the Cooper-pair binding energy as 
function of — |J7| for different values of the band filling n. 
A typical result for n — 0.4 is shown in Fig. |^. The results 
show that the binding energy is an increasing function of 
\U\ both for q — I (corresponds to zero field) and for 
q = 2,3. Neglecting spin-splitting, the binding energy is 
clearly higher for g > 1 and is an increasing function of q 
for the extreme fields we are considering. An increase in 
q makes the allowed energy bands narrower, and thereby 
increases the density of states, and hence this result is 
to be expected. Since the Fermi level normally lies in a 
band, the density of states at the Fermi level is higher 
than in zero field. 
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FIG. 1: Binding energy for Cooper pair as function of on- 
site effective potential for different magnetic fields (q). The 
Zeeman spin splitting is neglected and the band filling is n = 
0.4 (0.5 is half filled band). 

When one includes the Zeeman spin splitting, the Hof- 
stadter bands will split. As the effective Lande g- factor 
increases, the spin up and spin down bands will move 
relative to each other, and ultimately be completely sep- 
arated. A typical result for U = —At is shown in Fig. |^. 
The binding energy is not a monotonically decreasing 
function, but increases and decreases as the bands over- 
lap more and less. When the bands have passed each 
other and the spin up states are increasingly separated 
from the spin down states the binding energy is a de- 
creasing function of the factor. 

A similar plot as in Fig. |] is shown in Fig. ^. Here the 
binding energy as function of — is plotted for different 
values of the spin splitting for g — 3. 
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FIG. 2: Binding energy for Cooper pair as function of the 
efi^ective Lande p-factor (spin-splitting). The band filling is 
n = 0.4 and g = 3. 
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FIG. 3: Binding energy for Cooper pair as function of on-site 
effective potential for different values of the effective Lande 
g-factor for band filling n = 0.4 and g = 3. 



before the decrease towards the normal Hc2 from above. 
If the effective Lande g-factor is not too high, this result 
is still valid when we include the Zeeman effect. The 
results tell us that in the model we have considered, 
re-entrant superconductivity as a function of magnetic 
field will appear if the on-site attractive potential is suf- 
ficiently strong. For real materials, phase transitions to 
different phases may appear for lower fields than those 
we have studied. 



IV. CONCLUSION 



For the attractive Hubbard model in extreme magnetic 
fields the above results indicate that the binding energy 
is an increasing function of g, i.e the binding energy is 
higher in a magnetic field (extremely strong) than in the 
zero field situation. It is not clear what will happen for 
even higher values of g, but there has to be a maximum 
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APPENDIX A: TRANSFORMATION OF THE INTERACTION TERMS 

By using Eq. ^ and introducing rj = Vi + 6, we can write 



2 

:.],rT.a' 



1 ^ ^ y^|^|^gj(K4-K3)-5gi[(K-2+K4)-(Kl+K3)]-ri 
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X 



Et t 



Since Uk,,i{x) is periodic with period q in the x-direction (and 1 in the y-direction), we can use a Fourier-transform 
over the magnetic Brillouin zone 

= ^$]<KG.)e^^^^% (A2) 

or 

u^Ar^) = ^ E ^-.^ (G)e*^-'-' , (A3) 

where Gx — 0, 27r/g, •••((/— l)27r/g and Gy = 0. By using this, we get 

^^^^ E E ":.A(Gi)Sf..A(G2)xS:3.,3(G3)u.,,,,(G4)Eni^l)e^<'^^-'^^)-'e^(^-«3)-^ 

2 



= JVl5(„2+Gf2) + («4 + G4).(~l+G;i) + (»3 + G3) 

If we then define Kq by 

(/«i + Gi) - (k2 + Ga) = (k4 + G4) - (k3 + G3) = /tg, (A5) 

and 

v{Kq)^Nj2ymy'^'''^ (A6) 

we can write 

^^"2^ E ^('^9) E "iJ^2+f^,-(Gi-G2),£i(^l)^'^2.«2(G2)u^,3/3(G3)u^3+«,_(G4-G3),<!4(G4) 



X 



E '^l2+K,-(Gl-G2)/l,0-'^'«2/2,0-ci3,^3,cr'CK,3 + K,-(G4-G3),^4,0-'- (A''') 

This can be written on a more symmetric form by letting K2 k- — ^ , k' — 

^^ = ^ E ^('"9) E ^;^+^-(G4-G2)A(^l)''--^.^2(^2)«;,_:^,3(G3)iI,, + ^_(C^_G3)A(<^4) 
' Gi---G4 

E'^L+^-(Gi-G2)A,<T'^'^-^/2,^^l'-^,f3,^'^«' + ¥-(G4-G3)A,'T'- 
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Now we perform a BCS like truncation 



cr = —a 

k! — —k; (— Tr/q, — tt) < k < (7r/g,7r), 



(A9) 



where n is in the magnetic Brillouin zone, i.e. we assume that the Cooper pair is a spin singlet. Due to the 
antisymmetry of the spin part of the wave-function, the space part has to be symmetric. If /«' + /t 7^ 0, the wave- 
function for the pairs center of mass would have a modulation. This modulation will generally be incommensurable 
with the underlying lattice. Thus 
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By letting k + — > k, and k -t- ^ ^ we get 

= E E ^;^-(Gi-G2)a('^i)^'^'.^2(G2)w1«^^3(G'3)u_^'-(G,-G3),£,(G4) 

cL(Gi-G2)A,'tC'''A,<tcL^^£3^_<,C_«/_(G4-G3),^4,-<T- (AH) 



Gi---G4 
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To transform this back to the magnetic Brillouin zone we use Eq. 12c, in addition to the fact that the combination 
g-iK-T-iy* ^(xi)(^^ ^ ^ has to be invariant if we let k ^ /« + G (as can be seen from Eq. |^). Then 



We can then write Eq. All 



as 



= ^ E ^ ~2 E <A(''l)"'^'A(''2)wl^,£3(r3)W-K',£4('^4) 



(A12) 



(A13) 



K,,K, ,cr 



Gi ■■■G4 
ri ■ ■ ■r4 

^ g-(Gi-G2)-rig(G4-G3)-r4gGi-ri-G2-r2+G3-r3-G4-r4 

^ ''K,^1,ct'''^'/2,0-cL^ £3 _g.C_K',£4,-0■ 



(A14) 
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X „G2-(ri-r2) G3-(r3-r4) t r , e C , c 

^ ^K/l,cr^K' ^^2,o■^-K/3,-o•^-K',f4,-cr 

i V{k - k') ulj,^ {r)u^'j^ (^)"l«/3 A (r')cl jj,^c^',f2,<Tcl^,^3,„^c_ 

K,k'.(T T-J-' 
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(A15) 



G.G' 



From Ref. ^, we know that the term V{k — k') can be written as 



(A16) 



where A = (t/, V, W, V, W), and 



B2{k 

Bi{K 
B^iK 



= Vn 

= Vn [cos(kk) + cos(Ky)] 
= V2N cos{kx) cos^Ky) 

= VN [cOs(k2:) ~ COs(Kj^)] 

= V2N sin(K2;) sin(Kj,). 



(A17) 



The upper Umit ry < 5 in the sum is due to the even symmetry of the space wave-function. We can then write 
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(A18) 



=B 



ri.G,G' 



To simphfy notation it is useful to introduce the shorthand notations ^ = G, G'}, and x = {^ii ^2, 1^} which allows 
us to write 



(A19) 



where = A,,. 
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